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Abstract 

Products of random 2x2 matrices exhibit Gaussian fluctuations around almost surely 
convergent Lyapunov exponents. In this paper, the distribution of the random matrices is 
supported by a small neighborhood of order A > 0 of the identity matrix. The Lyapunov 
exponent and the variance of the Gaussian fluctuations are calculated perturbatively in A 
and this requires a detailed analysis of the associated random dynamical system on the 
unit circle and its invariant measure. The result applies to anomalies and band edges of 
one-dimensional random Schrodinger operators. 


1 Main results 

The analysis of random finite difference equations such as harmonic chains or the one-dimensional 
Anderson model naturally leads to study products of random 2x2 matrices. Asymptotically 
these products converge to a Gaussian process which is characterized by the Lyapunov exponent 
and its variance EQ. Often one is interested in a perturbative regime of small randomness. For 
example, the zero frequency limit of a harmonic chain is of this type, as well as the analysis of 
band center and band edges of the Anderson model (Ishl IKW. IDG] . Also the behavior near the 
so-called critical energies of a random Kronig-Penney model is of this perturbative nature jDKSj . 
It is the object of this work to develop a rigorously controlled pertubation theory for both the 
Lyapunov exponnet and its variance within a generic model covering all the situations alluded 
to above, and potentially others. 

More precisely, let us consider a family (7A, n )n>i of random matrices in S1(2,R) which is of 
the form 

Tx,n = exp(A'P n + A 2 Qx,n) , (1) 

where A > 0 is a small coupling parameter and V n , Qx, n £ sl(2, M). The matrix Qx )U is supposed 
to be analytic in A, and V n and Qx, n are independent and identically distributed on a bounded 
set. The distribution of 7a, n is thus supported in a neighborhood of size A around the identity. 
Averaging over this measure is denoted by E. Associated to any random sequence of matrices 
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are products J^Li 7a , n = 7a, jv • • • 7a, 1 and their assympototics satisfy a 0-1 law characterized by 
the Lyapunov exponent 7 a (EU It can be calculated by 


7a 


N 


lim — E 

N—too N 


logdi 


A,n^| 


n=1 


( 2 ) 


with some initial condition e G M 2 , ||e|| = 1. Furthermore, a central limit theorem [iTutl IBL] 
states that the expression 

■^= ^ lo g(ll n^ell) - ^7 a^ (3) 

converges to a centered Gaussian law with a variance denoted by G\. Both the variance and the 
Lyapunov exponent are independent of the initial condition e. Clearly the distribution of V n 
should have a crucial influence on the values of 7 a and a\. The following theorem summarizes 
the main results of the paper. 


Theorem 1 Let the i.i.d. random matrices 7a ,n be of the form (JT]) with the distribution supported 
on a bounded set. Furthermore, it is supposed that in each of the cases below a certain linear 
combination of the entries of V n specified below has a strictly positive variance. 


(i) If the averaged perturbation is elliptic in the sense that det(E[7 ? ri ]) > 0, then there is 
a constant C e , which is calculated explicitly from the variances of the entries of V n in 
Sectional such that 

7 a = C e A 2 + 0(A 3 ) , cr A = G e A 2 + (7(A 3 ) . 

(ii) If the averaged perturbation is hyperbolic in the sense that det(E[7 > „]) < 0, then with C h = 
a/— det(E(7 > n )), 

7a = C h \ + 0( A§) , ua = 0( A§) . 

(iii) If the perturbation is centered EJTA] = 0, then there exists constants C s and C' s such that 

7a = C S A 2 + 0(A 3 ) , cr A = a B A 2 + 0( A 3 ) . 


There have been numerous works in both the physics and the mathematics literature on 
products of random matrices of the form (JTJ). One is an influential work of Kappus and Wegner 
[KW ] on so-called band center anomalies for the Lyapunov exponent at the band center of the 
one-dimensional Anderson model. The term anomaly reflects that the standard perturbation 
theory breaks down and that this leads to an enhancement of the Lyapunov exponent. These 
anomalies and also the band edge were then further studied by Derrida and Gardner DS3- 
In the sequel, many rigorous works analyzed the Lyapunov exponent of the Anderson model 
in these particular citations [BK1ICK1ISVW1ISB1 155] . In [SS] any model of the form ([!]) was 
called an anomaly and such anomalies were further distinguished to be of first order elliptic if 
de^EfTA]) > 0, of first order hyperbolic if det(IEfT^rJ) < 0, and of second order if E[7 \\ = 0. 


2 






















This terminology will be used below and corresponds to the three cases in the theorem. Let us 
point out that the parabolic first order case det(E[7 :> rt ]) = 0 and E [P n ] ^ 0 is not covered by the 
theorem, but can be considered as non-generic. 

The main novelty of Theorem [1] are the formulas for the variances. In fact, the Lyapunov 
exponents in (i) and (iii) were calculated in [SS] and m respectively, and the case of (ii) was 
sketched in [ SSl DKS], but there are actually considerable technical difficulties to make the 
analysis rigorous, see Section [H On the other hand, no other work on variances at anomalies is 
known to us. The paper |SSS] developed a perturbation theory for the variance in situations where 
the lowest order matrix 7o, n is a non-trivial rotation rather than the identity. This corresponds 
to energies away from the band center and the band edges in the Anderson model. To control 
the variances in Theorem [1] it is necessary to go far beyond the analysis of [ SSS j. 

In the first order elliptic case of Theorem [Dp), the result implies the equality <j\ = 7 ^ + C)(A 3 ). 
Hence the asymptotic distribution of the product of random matrices is described by a single 
parameter, up to errors of higher order. In the framework of random Schrodinger operators this 
is referred to as single parameter scaling. It was already shown to hold away from anomalies in 
[SSSj . but only to lowest order in perturbation theory. This is now extended to the first order 
elliptic regime. For an Anderson model, this covers energies inside of the band, but close to 
a band edge (see Section [6]). On the other hand, energies outside of the band, but close to a 
band edge, correspond to a first order hyperbolic anomaly as in Theorem [TJii). In this situation, 
the Lyapunov exponent is given by its deterministic value, up to fluctuations which are of much 
smaller order as expected (because the system is a random perturbation around a hyperbolic 
one, albeit a very weakly hyperbolic one). In particular, there is no single parameter scaling in 
this regime. We are not able to determine whether <j\ is of even lower order than 0( As). In the 
centered case of Theorem [D(iii), the Lyapunov exponent and variance are of the same order of 
magnitude, and the formulas in Section [5] show that an equality C s = C' s is to be considered a 
coincidence so that single parameter scaling does not hold in this case. 

In the final Section [6] of the paper, it is shown how Theorem |T] can be applied to the Anderson 
model as well as two other models, namely a harmonic chain and the Kronig-Penny model. Both 
have already been considered in the work of Ishii |lshj . but there has been continuous and also 
recent interest in them ( e.g. EH! and [ DKS ] respectively). The remainder of the article is 
mainly devoted to proofs. Apart from perturbative formulas and recursive arguments, the main 
technical tool is the analysis of singular differential operators. 

Acknowledgements: We thank the referee for careful reading and several comments that 
improved the presentation. This work profited from financial support by the DFG. 
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2 From random products to a random dynamical system 

This section contains some preparatory material, most of which is already contained in |BL . SSl 
ESS], but is needed to fix notations and make this work self-contained. 


2.1 Dynamics on the unit circle 

For any P = (°_ b a ) € sl(2, R) one has P 2 = — det(P)l 2 which implies 


. . , ., ,, sinh(Ad) 

exp(AP) = cosh (Ad) 1 2 + -——-P 

C i 


d = \J — det(P) e C . 


A dynamics 6 e S P (9) e on §* = [0, n) is induced by 


e s XP {B) ~ 


exp(A P)ee 
|exp(AP)e e || 


ee = 


cos(6>) 

sin(d) 


This dynamics can readily be calculated in terms of harmonics in 9, namely with v = (j) 


0 2iS\p(j)) _ 


exp(XP)e g cosh (A d)e ie + sinh ] Aa!) v *Pe 6 


v*exp(\P)e e cosh(Ad)e - * 6 ' + sinh ^ d) v*Pe e ’ 

so that in terms of the entries of P 

e 2 x{s XP {0)-0) = cosh(Ad) + s -^f^(ae- 2ld + f(c - b) + f(c + b)e ~ 2ie ) 
cosh(Ad) + sm M Ad ) ^ ae 2te — — £>) — ^(c + b)e 2ld ) 

Expanding in A therefore leads to 

Sxp(6 ) = 6 + A p(9) + 1 \ 2 p{6)d e p(6 ) + (9(A 3 ) , 

where, still for P = ( ° 

V c —a 

p{9) = —asin(26 l ) + |(c— b) + |(c + b) 008(26*) = —asin(29) — b sin 2 (9 ) + c cos 2 (6 ) ) . (4) 


All the above applies directly to 7a, n- However, in order to deal with the various types of 
anomalies, it will be advantageous to make an appropriate basis change M G Gl(2, R) from 7a , n 
to 

7a, n = MTx^M- 1 = exp(A P n + A 2 Qx,n) , 

with P n = MVnM- 1 and Qx, n = MQx^M^ 1 respectively. The adequate choice of M will be 
made for each type of anomaly further below. Let us denote the dynamics associated with 7 a, n 
by Sx, n = Sx,p n +xQ Kn - According to the above, 

SxA°) = 0 + A p n (9) + A 2 q n {9) + ^A 2 p n (e)d e p n (d) + (7(A 3 ) , (5) 
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where p n and q n are defined by ([4j) from the entries of P n and Q o, n respectively. Now given 
any sequence T\^ n or equivalently ( P n , Q\,n)n>i and any initial condition # 0 G , one obtains a 
sequence 

On+l = S X ,n(0n) ■ 

If P n and Qo,n are random, this provides a random dynamical system on the circle S^.. 

2.2 Invariant measure 

There is an invariant measure u\ on corresponding to the action of the family (TyrJneN, 
which is given by 


r m "x(dO) = E r f(S x , n (9)) Ux(d0) , / e C($l ). 

Jo Jo 

From the Furstenberg theorem [BL j follows that the invariant measure is unique and Holder 
continous, provided the Lyapunov exponent of the associated system is positive. This actually 
holds in all cases considered here. Of interest will be the following ergodic sums 

N—l 

h,M) = x eE-WG 

n =0 


and especially in their limit 

d\(f) = Jim )■ 

N-> oo 

This limit can also be expressed in the terms of the invariant measure v x for which several 
notations are used: 

W) = rm«W = n(f) = E „[/]. 

As in [SSSj the calculation of the variance a\ is based on a control of the correlation sum 


Uf) 


Ei 


n =1 


Here J\(f) depends on the initial condition 9 0 and E n denotes the expectation over all T\ >rn with 
m > n. Hence Ei = E. It is shown in [BL ] that positivity of the Lyapunov exponent implies that 
) is finite for Holder continuous functions. Let us point out that / (->• J A (/) is linear and that 
constant functions are in the kernel of this map. A quantitative bound is given in Proposition [2] 
below. 


2.3 The Lyapunov exponent and the variance 

Let us introduce a family (g n ) n >\ of random variables by 

9 n = ^logdlT^.JI 2 ). 
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The Lyapunov exponent can then be defined as 


N 


7a = 


lim T7 E V] 9n = I\(g) , 

TV—>-00 J\l * ^ 


( 6 ) 


77.—1 


where g : —> K. is given by g(9 n -i) = E n [g n ], while the variance is 


N 


ax = i^loN S (^- 7 a )(^- 7 a ) 


1 


n,k=l 

N 


lim — > E^. 

TV-^oo N ^ X 
n =1 


N—n 




9n9n-\-m Ta 


777— 1 


The positivity of the Lyapunov exponent implies again that the sum over m is convergent 
even if N —> oo. Furthermore, E n (g n+m ) converges exponentially fast to 7 a as m —* oo [BL'. ISSS'J 
and the summands in the sum over n converge in expectation so that (with integration w.r.t. v x 
over #o) 


o\ = E„ E 


9 2 i ~ 7a + 2 ( 9i9m ~ 7a) 

m =2 


= E„ A E 


gl - 7a + 2 91 E 2 ( 9m - 7 a )) 

m =2 


(7) 


= E„ A E [gl(d 0 ) - j 2 x + 2g 1 (e Q )Jx(g){e 1 )\ = E Ux [g 2 - ^ + 2E( 9l J x (g) o S^)] . 


2.4 Basic perturbative formulas and estimates 

Replacing the expansion 

\\T x , n egl = 1 + 2A e* e P n ee + \~{2e* e Q n ee + e* 9 P*P n ee — det(P n )) + 0(A 3 ), 
into the definition of g n leads to 

9n = Ae; n _ 1 P n e e „_ 1 + A 2 (e; n _ l Q n e 0n _ 1 + ^ det(P n ) - (e^Pneg^J 2 ), 

up to terms of order O (A 3 ). Using 

Tr(4) a-d 


e* e Ae e = 


a - d b + c . . (a b\ 

+ — cos(29) + — sm(29), A = ( c , 


one can rewrite g n in terms of the entires of P n = p 1,n ^ 2,11 ) and Q n = ^ ^ 1,n " 


P3,n Pl,r 


93,n 9l,r 


9n =A^Pi,„COs(20 n _i) + i(p 2 ,n + p 3 ,„) sin(20„_i)^ + A 2 ^i[4p 2 n + (p 2 ,n + P3,n) 2 ] 

+ [gi + ^(>3,n - Pgr,)] cos(20 n _!) + ^ [g 2 ,n + 93,n + Pl,n(P2,n ~ P3,n)} Sm(20 n _!) (8) 

- ^[4- (P2,n + P3,n) 2 ] COs(46»„_i) - ^[pi,„(p 2 ,n + P3,n)] sin(40 n _i) ) + C>(A 3 ) . 
Using these formulas, the following quantitative a priori estimate on J x (f) is proved in 
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Proposition 2 There is a constant C such that for any Holder continous function f with Holder 
norm ||/|| a; 


Uf) < 


c 


A7a 


3 Elliptic first order anomaly 

This section considers the elliptic first order anomaly and proves item (i) of Theorem CD As 
det(E["P n ]) > 0, there is a basis transformation by M such that the matrix E [P n ] is the generator 
of a rotation: 

p - ( 0 -V\ , p _ ( Pl,n P2,n~h 
71 U oi + " \P3,U + V -Pun 


(9) 


where rj = sj det E[P n ] and P n is centered. Thus to lowest order the dynamics is a deterministic 
rotation perturbed by random (centered) fluctuations: 

S\,n{0) = d + Xr] + X Q(p 2 ,n - P3,n) + ^(p 2 ,n + P3,n) cos(2 6) - p 1>n sin(26>)^ + 0(A 2 ). (10) 


Proposition 3 |SS] Up to corrections, the invariant measure at an elliptic first order anomaly 
is given for f G C' 2 (§J.) by the Lebesgue measure: 


w) = - / m<w + o(A). 


7r 


Proof. Expanding / G C ,2 (§) r ) in Taylor series and taking the expectation provides 

E„+i [f(S x , n (9 n ))] = f(6 n ) + X V f'(9 n ) + 0(X 2 ). 

By summing this up one obtains 

IxAf) = hAf) + A vhAf') + 0{X 2 ,N~ l ) , 

and since / — y fj f(9)d9 has an antiderivative, the statement follows. □ 

In particular, I\( cos(2-)) = J A (sin(2-)) = 0(A). This allows to calculate the Lyapunov expo¬ 
nent. Using (jHI) one gets from (jHJ) 

9n = A fpi, n COs(26> n _i) + ~(P2,n + P3,n) Sin(20 n _i)^ + A 2 (^[4^1 >n + (p 2 , n + P3,n) 2 ] 

+ [Ql + \{pi,n - pi,n) - \v(P2,n + P3,n)} COs(20 re _i) 

+ + q3,n + Pl,n(P2,n ~ P3,n ~ 2?7)] SUl(26> n _i) 

~\[^Pl,n - {P2,n + P3,n) 2 ] COs(46> re _i) - ^\pi, n (p2,n + P3,n)} Sin(40 n _i)) + 0(A 3 ) , 
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so that the Lyapunov exponent is 

7 a = 4" (4Var(pi) + Var(p 2 + p 3 )) + C>(A 3 ). (11) 

O 

Here we omitted the indexation with n, since the P n ’s are distributed identically. This formula 
also already appeared in [p5j . The calculation of the variance uses the following lemma. 

Lemma 4 Let f G C 3 + "(S) r ) with a > 0. Choose c G R such that the function f — c has an 
antiderivative F G C' 4+a (§).). Then, independent of the choice of F, 

Mf) = -^(f(»o) - MF)) + o(i). 

Proof. Let us begin by noting that Proposition [2] combined with y A = 0( A 2 ) as given by (TTT1) 
shows that J\(f) = 0( A -3 ) because / is a Holder continuous function. Furthermore, one has 
J\(f — c) = J\(f) so that it is possible to suppose c = 0 and f(9) dd = 0 from now on. Next 
let us recall the Taylor formula 

F(9 + e) = F(9) + e F'(9) + ^e 2 F"(9 + +) 

where 0 < e/e < 1. Replacing 9 = 9 n and e = 9 n+l — 9 n = S n ^\(9 n ) — 9 n and taking the average 
of (HOD leads to 

E n F{9 n+1 ) = E n F(S x , n (9 n )) = F(9 n ) + A r)F'(9 n ) + A 2 r{9 n ) , 

where r is some Holder continuous function because / G C l+a (§*). Replacing this into J\ and 
separating the first summand, one has 

Jx(F) = (. F(9 0 ) - E„ a (F)) + J X (F) + Xr)Jx(F’) + A 2 J A (r) . 

Therefore, for / = F ', 

Mf) = - T (F(e 0) - MF) + A 2 J A (r)) ( 12 ) 

As J\(r) = (9(A~ 3 ) by the a priori estimate, this implies J\(f) = 0( A -2 ). Up to now only 
/ <7 U 1 + "(S) r ) was used. As / G C 2+ "(§ 4 ), one has r G C 1 + “(S) r ) and hence also J\(r) = 0( A' 2 ). 
This implies J A (/) = C^A -2 ) for / G C 2+ "(§ 4 ). One further iteration is needed to calculate the 
value of J\(f). As / G C 3 + “(§A), one has r G C 2 + "(S) r ) so that J\(r) = 0( A -1 ). Replacing this 
in (1T2T) concludes the proof. □ 

Lemma H implies J A (cos(j-)) = --^sm(j9 0 ) + 0(1) and J A (sin(j-)) = cos (j9 0 ) + 0(1). 

For the calculation of the variance the index n is again suppressed: 

J\(g) = ^(-E[gi + i(^-^)]sin(20 o ) + ^E[g 2 + g3+Pi(P2-P3)]cos(20o) 

+ ^E[4p| - (p 2 +p 3 ) 2 } sin(46> 0 ) - ^E[pi(p 2 + p 3 )\ cos(46> 0 )) + 0( A 2 ) 
Replacing into dTJ) shows that 07 = y A + (9(A 3 ) with y A given by (HID . 
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4 Hyperbolic first order anomaly 

A hyperbolic first order anomaly corresponds to item (ii) of Theorem [D Now it is possible to 
choose a basis change M such that 


fPri. 1 Pn ,2 \ 
\Pn, 3 ~Pn, 1) 


E(p n , 2 ) = E(p W)3 ) = 0 , r] = E (p nA ) > 0 . 


(13) 


Then the dynamics is to lowest order (higher orders are irrelevant and thus suppressed in the 
following) 


S\,n(0) = 0 + X (~p n ,i sin(26>) - p n>2 sin 2 (6>) + p n;3 cos 2 (6>)) + 0( A 2 ). (14) 

Hence the averaged dynamics is simply E S\ >n (9) = 9 — X E(p n l ) sin(2$) + 0( A 2 ) which has 6 = 0 
as a stable and 6 = | as an unstable fixed point, and no other fixed points. Hence, under the 
averaged dynamics any initial point 9 0 ultimately reaches the stable fixed point 0, except if 9 0 = | 
is the unstable fixed point. The first aim is to show that this behavior is stable in the sense that 
the invariant measure v\ converges to a Dirac peak on the stable fixed point as A —» 0. This 
is actually a delicate endeavor as will become apparent shortly. Once u\ is analyzed, Lyapunov 
exponent and variance can be deduced in Section 14.21 

4.1 Random dynamics on S 1 with two fixed points 

Let us consider for a function F £ C 2 (§ 2 ) the Taylor expansion 

F(S x , n (6)) = F{9) + A p n (0)F\9) + 0(X 2 \\F \\ c .). 
where p n (9) = —p n ,isin(20) — p U) 2 sin 2 (6 l ) +p n ,3 cos 2 (9). Thus the associated ergodic sums satisfy 

N-l 

I\,n(F) = -£>[F(0 n )] = I X , N (F) + A/ A , iV (E[p n ]F / ) + 0(X 2 \\F\\ C 2, N -1 ), 

n =0 

so that for every function / = A4F in the range of the first order differential operator A4 : 
C 1 ^) C(§i), defined by 

(. MF)(6) = 7} sin(20) F'(9) , (15) 

one has /a,v(/) = 0 (A||F||c2, A^ _1 ). Thus, upon integration w.r.t. the invariant probability 
measure v\, one has v\{f) = (9(A||F|| C 2). Every function in Ran(A4) has zeroes of order at least 
1 in both fixed points, but Ran(A4) contains sufficiently many functions to conclude that the 
weight of v\ is of order 0( As) outside of A 4 -neighborhoods of the both fixed points. Indeed, let 
us choose F(9) = — cos(29)/2r/ so that f(9) = ( MF){6 ) = sin 2 (20). Thus J A (sin 2 (2-)) = (9(A) 
so that, using the indicator function xa on A = {9 £ S^. | sin 2 (26*) > A 2 } which is bounded by 
Xa(9) < A _ 5 sin 2 (20), one deduces v\(A) = O(X^). 

Remark Of course, it is possible to implement this line of argument for many other functions. 
For example, taking F{6) = sin(26 l ) one obtains / A (sin(4-)) = (9(A). Analogously, with F(9) = 


9 




sin((2 + 4 k)9) for k £ N, one obtains inductively I x (sm(4k-)) = 0( A). Similarly, one obtains 
I x (cos(4k-)) = 1 + 0( A) for fceN. Furthermore, using directly the dynamics (TT4|) 


1 N—l 1 AT—2 1 N-2 

nT, e6 " = jv E E 9 ”- A, 'v^ Esin(2< '" ) + °( A2 ' JV_I) . < 16 ) 

n=0 n=0 n=0 


so that /\(sin(2-)) = 0( A) and, by the above procedure, /a(sui(2&;-)) = 0(A) for all k £ N. o 

However, the above strategy does not allow to say what the distribution of ^-weight of 
balls centered at the two fixed points 0 and | actually is. For example, one cannot evaluate 
7a(cos(2-)) by this method (its value is then equal to Ja(cos((2 + 4k)-)) for k £ N) because the 
function 9 y cos(26*) takes different values at the two fixed points. To go further it is necessary 
to zoom the dynamics into the unstable fixed point | and show that the fluctuations drive the 
random dynamics out of the neighborhood of this point. This will be done by a change of 
variables induced by a special Mobius transformation, in a manner explained in [SS] : 



!) 


exp ( A 





!) 



0 P2 
0 0 


+ A 


Pi 

0 
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where the abbreviation p n j = pj was used. The dynamics induced by the r.h.s. will be denoted 
by S\,n ■ ;> Sj. and the associated random dynamics by (9 n ) n > o- All objects in this rescaled 

representation will carry a hat in the sequel. One finds from (l5jl 

S\,n(9) = 6 ~ ^Pn ,2 sin 2 (0) + A (^Pn ,2 sin 2 (9)-Pn,i) sin(20) + O(X^) . (17) 

Now the lowest order term (in A^) is centered, and the first term with non-vanishing mean is 
precisely of the order (here A) given by the square of the lowest order term. Hence one is a 
situation where the techniques of [SB], in particular Proposition 3(vii), can be applied. However, 
we present a self-contained argument here in order to show that almost all the ^-weight lies in 
a neighborhood of the stable fixed point 0. 

Proposition 5 Suppose E[p|] > 0. Let \b be the indicator function on B = {9 e | |0| > A^}. 
Then 

o x (B) = 0{ A3) . 


Proof: Underlying the above basis change to the rescaled dynamics (1171) is the diffeomorphism 
t x :§l^ Si defined by 

t x {9) = tan -1 (A - 2 tan(0)) , 

It zooms into the unstable fixed point | and keeps 0 conserved. Then S x>n = tf 1 o S X n o t x and 
one is led to study the random rescaled dynamics (6 l n )„>o given by 9 n = tf 1 (9 n ). The aim is to 
control the Birkhoff sums /at(/a) of the rescaled dynamics for the functions f x = / ° t x , which 
are readily connected to the original ones via 

N—l 1 N—l 

hAf) = t = hAf »)• 

n =0 n=0 
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With the aim of controlling such Birkhoff sums let us again expand for a given function F G 
C 2 (§^) using (ED and E (p 2 ) = 0: 

E F(SxM) = HO) + X^^{sm 4 (9)F ,, {9) + (sm 2 (9)-2k)sm(29)F l (9)^ + 0{\l\\F\\ c2 ) , 

where k = E(pi)/E(p|) > 0. Defining a first order differential operator C : C 1 (S^ r ) — > £7(8*) by 
(Cg)(9) = sm 4 (9)g'(9) + (sin 2 (0) -2k) sin(29)g(9) , (18) 

one concludes that the Birkhoff sums satisfy 

I\,n(F) = Ix,n{F) + ,n(CF') + 0(\%\\F\\ C 2, N -1 ) ■ 

Therefore the function p = CF' satisfies 

IxAp) = 0(\^\\F\\ C 2,N- 1 )- (19) 


Let us choose in (TT9|) 

p{9) = exp(— k cot 2 (#)) , 

then solve p = CdgF for F and show that this solution satisfies H-FHc 2 < G uniformly in A (note 
that F depends on A via p). Hence one can first solve the first order equation CG = p for G. 
There is actually a two-parameter family of solutions G in the neighborhood of the singularity 0 
due to Proposition 3(vii) in |SS] (because the order of singularities and the signs of the coefficients 
are precisely as required there). The first parameter is fixed by requiring G to be continuous on all 
§1 (periodicity) and the second one by adding an adequate multiple of the (smooth non-negative) 
ground state G 0 satisfying CG 0 = 0, explicitly given by 


Go(6) 


exp(— k sin 2 (9)) 
sin 2 (0) 


This assures that G has vanishing integral and therefore has an antiderivative F. This explicit 
solution G = F' of CG = p is (as in [SSj) 


F\9) 


exp(— k sin 2 (0)) 
sin 2 (0) 




exp (A; sin 2 (s)) 
sin 2 (s) 



( 20 ) 


The function F' and its derivative F" are obviously regular on an interval I C §* not containing 
9 = 0 with bounds on / that are uniform in A (because p < 1 uniformly in A). More delicate is 
the behavior at 9 = 0 because the integrand becomes singular. However, by de l’Hopital’s rule: 


„ „ Ci + ft p(s) sin 2 (s) exp(/c sin 2 (s)) ds 

lim F\9) = lim —- J ° 2 \ K 

exp(fc sin 2 (9)) sin 2 (9) 


lim — 

0->o (1 


p(9) sin 2 (9) 
cos(2 9) — 2k) cot(0) 
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which vanishes due to the special form of p. Next calculating F" from the differential equation 
CF' = /, 


lim F"{9 ) = lim 

9^0 d->0 


p(9) - (sin 2 (9) -2 k) sin(2 9) F'{9) 
sin 4 (0) 


Since both denominator and numerator tend to 0 for 9 —> 0, de l’Hospital’s rule can be applied. 
However, using it directly is not helpful and it is better to first divide both by (sin 2 (0) — 


2 k) sin(20)G o (0). Define p{9) = 


m 


(sin 2 (0)—2 k) sin(20) 


and p{9) = 


sin 4 (6>) 


(sin 2 (0) —2 k) sin(2 0) 


and use that 


'Go(e)- 


GoV) pW 


Thus one obtains 


lim F"{9) = lim — ^ ^ „ = 0. 

e->o 1 + p'{9) 


One concludes that the C 1 -norm of F' is uniformly bounded in A. Let us note that jSSl Propo¬ 
sition 3] presents an iterative procedure showing that also further derivatives of F' exist and are 
bounded. 

Now let us consider the transformed p(9) = fot^ 1 (9) = exp(— k 5 cot 2 (6 1 )). Then /a,at(p) = 
In(p) = (P(A^||F|| c 2 )- Next set A' = {9 G S 4 | \9 — || < A^}. Then there is a constant C such 
that xa' < Cp and therefore v\(A') < C'X^. On the other hand, it was already shown above 
that the set A = {9 G | sin 2 (29) > Az} satishes v\(A) = 0( A 2 ). Combining these estimates 
concludes the proof. □ 


Corollary 6 Suppose E[p|] > 0. For every function f G C 1 (S 4 ) 


Mf) = m + wile 1 a*), 

while for f G C 2 (§ 4 ) 

v\(f) = /( 0 ) + 0(\\f\\c>\h- 

In particular, every sequence {v\ k )k> 1 of invariant measures with A k —» 0 converges weakly to a 
Dirac peak at the stable fixed, point 0. 


Proof: Let B = {9 G | \9\ > A«} and B c = \ B. Then by Proposition [5] 


Mf ) = 0(||/||ooA5) + [ f ( 9 ) V \( d 9 ) . 

JB 

Expanding f { 9 ) in a Taylor series of order 1 in A implies the Erst claim. For the second, let 
us replace / by f ( 9 ) = f ( 9 ) — |/'(0) sin(46 l ) so that /'(0) = 0. As z/ A (sin(4.)) = 0( A) by the 
argument preceeding (TT6]l . one has n A (/) = u A (/) + 0( A). Now a Taylor series of / to second 
order in A leads to the second claim. □ 
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4.2 Evaluation of the Lyapunov exponent and its variance 

Corollary [6] can be applied to the function f{9 ) = cos(20) and implies Ja(cos(2-)) = 1 + O(X^). 
Furthermore, it was already shown using (TT6|) that J\(sin(2-)) = 0( A). This allows us to calculate 
the Lyapunov exponent based on the expansion ([HD up to order A: 

7a = Xrj + 0( A*). 

This shows the first formula in Theorem [Tfii) under the condition E[p|] > 0 on the variance of 
the entries of P n . The second claim in Theorem [Tf ii) is an estimate on the variance a x . This will 
be based on (JTj). As it is already known that c/i, g and 7 a are of order 0( A), it is thus sufficient 
to show that J x (g) = 0( A^), which results from the following 


Proposition 7 Assuming 9 0 = 0. Then for any f e C' 3+Q (§) with a > 0 

Jx (/) = 0( A-i). 

Proof. Again let us use the rescaled dynamics defined in (IT7|) . Let C be deffiied by (IIS]) and 
set /a = / o t\ as above. Furthermore, let v x be the invariant measure of the rescaled dynamics. 
Then v x (/) = v x (/) and there is a connection between the correlation sums of the two processes: 


us) = Ei [y (/(«„> -un) 

n=1 


n =1 


Ufx)- 


For / e C 1+Q (§ 3 ) with / 


CF', a Taylor expansion with some Holder continuous f gives 


Ja(F) = F{6 0 ) - hx(F) + J\(F) - \J X (CF') + A§ J A (f), 


so that 

Jx(CF’) = ^(F%)-v x {F)) +A hj x {r). 

Because 7 a = 0{ A), the a priori estimate gives J x (f) = 0(\~ 2 ) and one hence obtains J X (CF') = 

3 

0( A - 2 ). As mentioned before, every Holder continuous function with zero at 0 lies in the range of 
£. Furthermore, J\ is invariant under the translation by a constant. Thus, J x (r) = J x {f — f(0)). 
By (H9|) the invariant measure z >a is concentrated at 0, so the difference F(9 0 ) — v x (F) is of order 
A 2 because 9 0 = 0. Recursively it follows that J x (f) = for every function / e C' 3 +a (S 3 ). 

Hence for a function f\ = / o tf 1 one has 

Ufx) = Uf) = 0( A-3) . (21) 


Next for a given smooth function / let us define 

= /W-/(O)-(/(f)-/(O))sin(^ 1 (0))- C Asin(2^ 1 (0)), 


with the constant defined by 


ca = 


1 f m-m-vw-fmMtfm 

n J s i sin(2 9) 


d9 < 00 
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This constant is uniformly bounded in A as can be seen by analyzing the contributions around 
the singularities. For example, for small a > 0 one has J“ a < C and using that also 

[ a sm (t?(6)) f a ( A tan 2 (6 l ) V 1 Of tan 2 (g) V 1 

]_ a sin(2 9) J_ a \1 + Atan 2 (6*)y sin(2 6) J_ a \1 + tan 2 (6*)y sin(2 9) 


is bounded, one sees that also the second summand has a uniformly bounded contribution. 
Furthermore /(0) = /(|) = 0 and 


[ si n(2t A 1 (g)) = f X-2 

isi sin(2$) i s i Asin 2 (6») + cos 2 (6>) 



(^ 1 )'(0)d0 


7T, 


so that also / gl s /^ dd = 0. By (TT5T) one concludes that / lies in the range of M.. By (J21}) one 

now concludes J\(f) = J\(f) + 0( A - ^) and it thus only remains to analyze the correlation sum 
of / = A 4F. Expanding F e ( 7 2 +a(§i) in Taylor series and using the operator A4 defined in 
(USD, one gets 

Jx{F) = F(0 o ) - vx(F) + J\(F) - A ij J\(MF) + A 2 J A (r), 


for some Holder continuous function r, since F G C' 2+ °(S). By Proposition [2] and 7 a = 0( A), r 
satisfies the a priori bound J\(r) = 0( A -2 ). Hence 

Jx(MF) = -^(F(6 0 ) - v x (F) + 0(1)) , 

so that Jx(M.F) = (P(A -1 ). By the above applied to r one has J\(r) = Jx(r) + 0(X~ s) = 0(X~ l ). 
Thus actually Jx(JAF) = 0(X~ 5). This finishes the proof. □ 


5 Second order anomalies 

At a second order anomaly, the term of order 0( A) in (J5]) is centered. It is then not necessary to 
carry out a basis change M, but we nevertheless set P n = V n and Q n (X) = Q n ( A). Let us begin 
by recalling some results from [ SB, SS] which are needed below. First of all, one is naturally led 
to study the second order differential operator on C 2 (§^) defined by 

£ = |E n [p 2 ]<^ + E n [q n + \p n p' n ]de , 

because for F G C 2 (§^) one then has 

E n F(S n ,x(0)) = F(0) + A \CF)(6) + 0(A 3 ) . 

From this one can again control Birkhoff sums in the range of C. For sake of simplicity (weaker 
hypothesis are possible), it will be assumed that E n [p 2 ] > 0. This assures that C is elliptic. Let 
also introduce the formal adjoint C* = <9g|E n [p 2 ] — 8qE n [q n + ^p n p' n ]. Then [SSI, Theorem 3] 
implies the following result. Here scalar product and orthogonal complement are taken in L 2 (§) r ). 
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Theorem 8 Under the above ellipticity hypothesis, one has dim(Ran(£)- L ) = 1 and there is a 
unique smooth and non-negative function p G Ker(£*) with unit integral. For every function 
f G C 1 (§ 3 ) one has 

IxAf) = ( P\f) + 0{\ (A N)- 1 ), 

where ( p\f) = / §1 f(9) p(6 ) dO. 

7T 

Using Theorem |8]and (JHJ) as well as (|8l) , it is straightforward to write out the expansion for 
the Lyapunov exponent given in Theorem |T|(iii) with rigorous control on the error terms. The 
formula for C s was already given in |SB ] and invokes p which, in general, cannot be caluclated 
explicitly. For the calculation of the variance, one needs again to control J\(f). 

Lemma 9 For a function f G C 5+ "(§ 3 ) with a > 0 the correlation sum. J\(f) is given by 

■h(f) = -L(F(9„) - v x (F) + 0(A- 1 ) , 

with a periodic function F = £ _1 (/ — (/|p)) defined uniquely up to a constant. 

Proof. Let us first recall [HE] that the range Ran(£) includes all smooth functions that are 
orthogonal to p in the L 2 -sense. Lienee F is well-defined. Furthermore, J\ is invariant under 
shifts by a constant. Therefore J\(f) = J\(CF). As in the proof of LemmaHa Taylor expansion 
of F G C 3+Q (§ 3 ) shows 

Jx(F) = F(e 0 )-u x (F) + Jx{F) + A 2 J x (jCF) + A 3 J A (r), 

where r is a ffblder-continuous residual function for which an priori bound J\(r) = 0( A~ 3 ) holds. 
From this one deduces J\{CF) = 0( A -2 ) and hence Jx(f) = 0( A” 2 ). As all this equally well 
applies to r, namely J\(r) = 0( A -2 ), the equation actually already implies the claim. □ 

The calculation of a\ will again be based on (J7J) . By (JHJ) one has g n = A h n + A 2 f n + (9(A 3 ) 
with certain trigonometric polynomials h n and f n . By hypothesis, one has E [h n \ = 0. ffence 
the function g = ~E n [g n \ has an expansion given by g — A 2 / + 0( A 3 ) with / = E [/„]. Let 
F = C~ 1 (f — ( p\f )). Cancelling out the A-factors, one then concludes from Lemma 0 

Jx(g) = vx{F) - F(9 0 ) + 0( A). 


As 7 a = C^A 2 ), one thus has 

&x = E^Ei [A 2 h\ + (A hi + A 2 /i) Jx(g)\ + 0( A 3 ) . 

Replacing J\(g) one sees that due to E[/ix] = 0 the variance satisfies cr A = 0( A 2 ). In principle the 
above formula gives all the lowest order contributions, provided the contribution of order 0( A) 
to Jx(g) can be calculated. This provides C' s and it appears to be a coincidence for it to be equal 
to C s . 
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6 Examples 

6.1 Harmonic chain 


A harmonic chain consists of sequence of masses m n > 0 with displacements u n from it’s equlib- 
rium position for which the stationary equation of motion at frequency c o is 

YYl n 0J 'U'n ^n+1 n T U n —\ , Tl £ E . 

The masses are supposed to be i.i.d. randomly distributed with a distribution of compact support 
bounded away from zero. This is a Jacobi matrix and the associated transfer matrix is given by 

T _ /2 — uj 2 m n -l\ 

In,OJ - ^ x () J ■ 

The transfer matrix is in Sl(2, M) and it is conjugate to a rotation. The aim will now be to 
study the Lyapunov exponent and variances of the associated random products perturbatively 
in oj. For small oj , T n ,iu lies in the vicinity of the Jordan block. Indeed, by the basis change with 

Mi = \ 0 ) one obtains 


MiTn^Mi 1 = 


1 l\ Ml 0 
0 f )- m ^ 1 0 


Now one blows up the vicinity of the stable point by conjugation with the matrix M 2 = 


oj 0 

0 1 


M 2 MiT n , u (M 2 Mi ) 1 — 1 + oj 


t 0 

1\ 

2 ( 

OJ 

Oj 

- oj z 

\-m n 

\ 


m n 0 

0 0 ;' 


The last basis change is done with the matrix M 3 = f \ZM^n\ 0 


T n ,u — MzM 2 MiT n , u (M 3 M 2 Mi) — 1 + oj | m n q 

\7 E [ ?i 


0 -^/E[m n ]\ 2 fm n 0 

m n n — OJ „ 


0 oy ' 


This is a first order elliptic order model already in the form of ((9]). Therefore the results of 
Section [3] giving the constant in Theorem []Ji) shows 

w 2 E [ml\ - E [m n \ 2 _ 3 _ 3 

7 " T-IR]- + ° {u ) = ff “ + >■ 

The formula for the Lyapunov exponent is known at least since the work of HE]. 
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6.2 Band edges of the Anderson model 

The one-dimensional discrete Anderson model is described by the stationary Schrodinger equation 


E'tpn 1pn+l T Ipn —1 T Xl> n ■ 


( 22 ) 


Here i/} n are the probability amplitude, E E R is the energy and the v n E K are i.i.d. potential 
energies according to a centered distribution of compact support. The transfer matrix is given 
by 


T = 

• n,cj 


E — \v r 

1 


The case E = 0 corresponds immediately to an anomaly once one takes two consecutive transfer 
matrices as the basic building blocks [KW . SB j. Here we will rather focus on the band edges 
| if | = 2, and show that this also leads to an anomaly [DGlISSj. Let us focus on a situation where 
E = 2 + w A for some w E R, again with the aim of calculating the Lyapunov exponent and 
variance perturbatively in A. Their behavior depends strongly on the sign of w. If w < 0, one 
has an elliptic first order anomaly. From the technical point of view, this problem can be dealt 
with the same procedures as the previous one substituting cu 2 by A and E[mJ by —w, and this 
leads to 


7J = AEMHA + 0(A !) = + 0( a!) 


w 


For w > 0 one has a hyperbolic first order anomaly. The basis changes into the normal form 


are consecutively given by the matrices M\ = 


0 


, M 2 = 


\2 0 
0 1 


and M :i — 


W2 

1 

— W 2 



respectively. Now the results of Section 0] show 


11 3 

7a = A2u?2 4- 0 (\*), 


= 0 (\i) . 


6.3 Random Kronig-Penney model 


The random Kronig-Penny model describes a quantum particle by a continuous one-dimensional 
Schrodinger operator with singular random potentials on Z. This model also reduces to the study 
of products of random transfer matrices, as described in detail in [DKSj . Let us just recall from 
I )1\SI the special from of the transfer matrices close to the so-called critical energies Ei = (wZ) 2 , 
l G N, of the system, namely for e > 0 and after adequate basis changes 


T^~ £ = R i 

n -rje 2 


1 + 


and 


2y/2vEi 


-1 1 

-1 1 




0 1 
1 0 


£ V. 


2 Ei 


n I 0 1 

1 0 


+ 0(e*), 


rpj^i 

^ n 


Ei+e 


1 — r ]£2 
0 


1 + 7/£i 


+ £2 


2y/2vEi 


+ s 


AE, 


+ 0{£2) 


Here the v n G M are again i.i.d. with a distribution that is assumed to have compact support 
(not concentrated in a single point) with expectation value v and deviations v n = v n — v. Finally 
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R i is a rotation with an angle 77 = y Now T El £ is again an elliptic first order anomaly 
with £5 playing the role of A. Then the Lyapunov exponent y El ~ e and its variance a El ~ £ are np 

3 

to the order £2 given by 


7 


Ei-e 


- v 2 ) 

16 vEi 


e + 0(e§) 


a El ~ e + 0{e i). 


On the other hand, T El+e is a hyperbolic hrst order anomaly and thus the above results imply 

J + O(J), a El+e = O(ei) . 
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2 £, 


The Lyapunov exponents were already obtained by the authors in the previous work [DKS ], but 
the proofs and the error estimates there are given too briefly. 
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